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摘 要 : 近年 来 ,灰色 系统 理论 在 各 个 领域 中 的 应 用 愈 发 广泛 .在 灰色 预测 模型 
中 ,GM(1,1) 模 型 最 为 核心 与 基础 ,但 该 模型 对 振荡 序列 难以 进行 模拟 , 且 振 荡 序 列 往 
往 无 法 通过 准 指 数 律 检 验 ,这 削弱 了 GM(1,1) 模 型 的 应 用 范围 .本 文 利用 GM(1,1) 模 
型 对 单调 序列 的 模拟 ,对 振荡 序列 的 1-AGO 序列 建立 GM(1,1) 模 型 .再 通过 引入 修 
正 项 的 IAGO 操作 ,还 原 到 对 原 振荡 序列 的 模拟 .建立 了 改进 的 GM(1,1) 模 型 ,弥补 了 
传统 GM(1,1) 模 型 对 振荡 序列 模拟 的 缺陷 . 
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= 1 引 


灰色 系统 理论 自 1982 年 由 邓 聚 龙 教授 提出 后 , 距 今 已 经 有 近 40 年 的 历史 了 , 其 起 源 于 相 
据 数据 量 较 小 的 信息 作出 预测 与 分 析 03. 由 于 灰色 系统 理论 集 时 间 响 应 模型 、 灰 色 关 联 分 析 
en 方法 和 灰 数 理论 于 一 身 , 能 够 对 实际 的 问题 进行 多 方面 、 多 角度 地 分 析 , 所 以 迅速 地 成 为 经 
LO 济 、 信 息 、 决 策 和 控制 领域 中 十 分 重要 的 研究 手段 . 其 中 ，GM(1,1) 模 型 具有 十 分 清晰 的 数学 
ee 图 像 ,良好 的 函数 分 析 性 质 ,以 及 它 体现 了 灰色 系统 理论 建 模 的 精髓 ,所 以 成 为 了 应 用 最 广泛 、 
CN 最 基础 的 模型 . 而 其 他 的 多 种 多 样 的 模型 , 如 GM(2,1) 模 型 、Verhulst 模型 等 等 , 都 是 要 么 通过 
CQ 蔡 换 GM(1,1) 模 型 中 的 灰 作用 量 得 到 #9， 要 么 将 GM(1,1) 模 型 的 时 间 响 应 函数 和 别 的 函数 看 
到 合 得 到 9) 下 面 就 GM(1,1) 模 型 ,进行 简要 的 介绍 . 
一 源 于 灰 导 数 的 定义 , 对 原 序列 x (进行 1-AGO 生成 , 得 到 x (有), 该 1.AGO 序列 的 灰 
‘© 导数 是 原 序列 .以 GM(1,1) 模 型 为 例 ,其 白化 方程 为 
[em dx (k 

= i 人 


az0(K) =b (1) 


其 中 zk) 二 0.5x0(k) 二 0.5xu(K 一 ID) 是 对 1-AGO 序列 进行 平滑 性 处 理 后 得 到 的 背 
景 值 序列 . 对 1-AGO 序列 建 模 , 则 (1) 式 等 号 左边 第 一 项 对 应 的 便 是 原 序列 x (Kk). 这 样 使 
得 方程 的 求解 变 得 方便 . 

由 (1) 式 解 出 的 时 间 响 应 序列 为 

xV (k++1) = [本 (GD 一 2 多 直 2 (2) 


其 中 的 参数 a、b， 可 以 采用 最 小 二 乘 进行 估计 ， 具体 步 又 可 参考 文献 00. 

但 在 实际 应 用 中 ,该 模型 对 原 序列 的 有 很 多 严格 的 要 求 ， 比 如 由 于 GM(1,1) 模 型 用 到 的 
时 间 响 应 函数 是 指数 型 函数 ， 所 以 为 了 保证 模型 的 准确 性 ,要求 建 模 的 序列 必须 满足 准 指 数 
律 . 然而 很 不 幸 ,， 现实 中 的 振荡 序列 , 一 般 都 难以 满足 准 指数 律 ， 从 而 导致 GM(1,1) 模 型 的 
相对 误差 很 大 . 并 且 , 指数 型 函数 构建 的 模型 始终 无 法 很 好 地 避免 单调 性 所 带 来 的 影响 ,， 特 
别 是 对 振荡 序列 的 模拟 结果 不 尽 人 意 ， 指 数 模型 在 这 方面 的 缺陷 已 经 有 许多 工作 证 明史. 
正 是 基于 这 种 情况 ， 本 文 提出 了 一 种 基于 高 阶 AGO 序列 建 模 的 办 法 ， 即 对 原 序列 进行 两 次 
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甚至 多 次 AGO 
让 数据 序列 更 力 


生成 之 后 , 再 进行 建 模 , 并 最 终 还 原 到 原 序列 的 方法 . 由 于 高 阶 AGO 生成 会 
0 趋 于 指数 分 布 , 从 而 使 得 GM(1,1) 模 型 的 有 效 性 得 到 很 大 的 提升 ,比如 我 们 


随便 取 一 个 序列 {50810,46110,51177,93775,110574,110524} ， 其 1-AGO 序列 为 
148097, 241872, 352446 , 462970} ， 绘 制图 像 作为 对 比 


{50810, 96920， 


图 1 某 数据 序列 与 其 1-AGO 对 比 图 


图 1 的 (a) 是 原 


序列 , (b) 是 其 1-AGO 序列 , 可 见 AGO 生成 的 确 能 让 序列 更 加 趋 近 于 指数 分 布 . 


以 及 , 高 阶 AGO 模型 在 模拟 还 原 时 引入 的 带 有 正 负 性 的 修正 项 , 一定 程度 上 弥补 了 指数 型 
函数 模型 建 模 机 制 上 的 缺陷 ， 即 单调 性 对 模拟 带 来 的 限制 . 


2 高 阶 AGO 生成 的 模拟 优势 


GM(1,1) 的 


时 间 响 应 序列 是 指数 型 的 ， 所 以 适合 描述 整体 呈现 单调 性 的 序列 ， 这 也 正 是 


其 难以 描述 振荡 序列 的 原因 . 并 且 振 荡 序 列 往往 无 法 通过 准 指数 律 检 验 . 考虑 到 AGO 生成 
能 使 序列 的 发 展 趋势 趋 于 指数 型 05，GM(1,1) 模 型 对 原 序列 建 模 效 果 较 差 ， 却 很 可 能 能 够 很 
好 地 描述 其 AGO 序列 . 


下 面 以 两 双 


昌 数 据 举 例 来 说 明 这 个 问题 . 取 振 荡 序列 {9,6,4,5,3} 和 {8,1,6,3,4)}, 对 


其 进行 GM(1,1) 建 模 以 及 对 其 1-AGO 序列 进行 GM(1,1) 建 模 ， 比 较 相 对 误差 . 


表 1 对 序列 建 模 与 对 其 1-AGO 序列 建 模 的 比较 


序列 模拟 值 相对 误差 1-AGO 序列 模拟 值 相对 误差 
9 \ \ 9 \ \ 
6 5.748612 4.19% 15 15.594331 3.96% 
4 4.814434 20.36% 19 18.850539 0.79% 
5 4.032065 19.36% 24 22.786666 5.06% 
3 3.376834 12.56% 27 27.544686 2.02% 
平均 相对 误差 14.12% 平均 相对 误差 2.96% 
序 负 模拟 值 相对 误差 1-AGO 序列 模拟 值 相对 误差 
8 \ \ 8 \ \ 
1 2.829706 182.97% 9 10.506067 16.73% 
6 3.247475 45.88% 15 13.530107 9.80% 
3 3.726922 24.23% 18 17.424579 3.20% 
4 4.277153 6.93% 22 22.440026 2.00% 
平均 相对 误差 65.00% 平均 相对 误差 7.93% 


由 于 出 于 对 比 效果 的 目的 , 这 里 暂 不 考虑 模型 前 检验 而 直接 建 模 (实际 上 它们 的 1-AGO 
序列 均 能 通过 GM(1,1) 建 模 前 检验 ) . 从 表 1 可 以 看 出 , 直接 对 振荡 序列 进行 GM(1,1) 建 模 
的 平均 相对 误差 非常 大 ,而 对 振荡 序列 的 1-AGO 序列 进行 GM(1,1) 建 模 ， 其 平均 相对 误差 基 
本 都 在 可 接受 范围 内 . 这 说 明 GM(1,1) 模 型 更 适合 对 振荡 序列 的 AGO 序列 进行 模拟 , 而 不 
是 直接 对 原 序列 进行 模拟 . 
由 于 GM(1,1) 模 型 本 身 就 会 对 准备 建 模 的 序列 做 1 次 AGO 生成 , 我 们 这 里 是 在 GM(1,1) 
之 前 就 对 数据 进行 了 1 次 AGO 生 成， 整个 流程 中 出 现 了 两 次 AGO 生成 . 我 们 将 这 种 做 AGO 
生成 次 数 超过 1 次 的 GM(1,1) 模 型 称 为 : 基于 高 阶 AGO 生成 的 GM(1,1) 模 型 , 在 下 文中 简称 
“高 阶 AGO 模型 ”. 

接 下 来 我 们 需要 解决 如 何 将 对 AGO 序列 的 模拟 还 原 到 对 振荡 序列 本 身 的 模拟 ,以 及 这 
种 模型 如 何 对 序列 做 出 预测 . 我 们 先 从 理论 上 分 析 该 模型 应 当 如 何 还 原 对 振荡 序列 的 模拟 . 


3 模拟 还 原 方法 


设 原 序列 为 x (Kk)， 其 预测 值 记 为 XO(Kk)。x 中 (有) 的 1-AGO 序列 为 x (Kk), 其 预测 值 
记 为 站 (Kk)， 以 此 类 推 . 原本 的 GM(1,1) 模 型 利用 IAGO 操作 
XV KF) 一 xO(FE 一 1 


二 Ee (i) 一 Sxo (0) =xO(D G3) 


实现 了 从 对 zx (6) 的 模拟 还 原 到 对 x@ (k) 的 模拟 . 但 由 于 建 模 之 后 ，(3) 中 的 1-AGO 序列 要 
取 模 拟 值 ， 这 便 是 GM(1,1) 的 误差 来 源 . 
如 果 我 们 对 1-AGO 序列 x (Kk) 进行 GM(1,1) 建 模 , 得 到 2-AGO 序列 x(k) 的 模型 . 想 
要 还 原 到 x (kK) 序列 就 需要 两 次 IAGO 操作 
GD x KD)— (x koD—x (一 2) 
一 xDO(F) 一 xD(KE 一 D =x (FR) 


(4 中 的 x 中 (如 序列 是 模拟 值 ， 而 x" (6 是 模拟 值 计 算出 来 的 . 这 时 再 对 x (有 ) 进行 计 

算 会 造成 误差 的 二 次 累积 ， 从 而 导致 模型 误差 可 能 比 直接 进行 GM(1,1) 建 模 更 大 . 所 以 需要 

一 种 办 法 来 降低 第 二 次 IAGO 操作 的 误差 .从 而 更 精确 地 还 原 到 对 原 序 列 的 模拟 ,让 高 阶 

AGO 模型 的 优势 得 以 发 挥 到 最 大 . 

设 序列 与 模拟 值 之 间 有 残 差 A 三 各 (4) 一 x (6)， 以 及 相对 误差 5 = |Al/x(k). 可 根 
据 A 的 取 值 情况 ， 取 模拟 修正 项 

xOR) = 60(k), A>0 

{ xzo( 昌 =G+50io(O，A<0 


(4) 


(5) 


如 果 取 消 相 对 误差 中 的 绝对 值 符号 , 那么 (5) 中 两 个 式 子 可 以 统一 为 


XVK) = I— 6) x kK), 6:=A/xd (Kk) (6) 


到 


然后 利用 x'@ (6 对 xz (Kk) 进行 计算 ， 从 而 将 对 x 中 (1) 的 模拟 还 原 到 对 x 中 (CK) 的 模拟 
{ xo(D xO KD 一 xO( 有 


ET WD 


我 们 仍然 以 振荡 序列 {9, 6,4, 5,3} 和 {8, 1,6,3,4} 作 为 例子 来 研究 这 个 方法 的 效果 . 


表 2 模拟 还 原 效果 示意 表 


序列 GM(1,1) 模 拟 ”相对 误差 高 阶 AGO 模型 模拟 值 相对 误差 
9 \ \ \ \ 
6 5.748612 4.19% 5.976452 0.39% 
4 4.814434 20.36% 4.022365 0.56% 
5 4.032065 19.36% 4.939851 1.20% 
3 3.376834 12.56% 3.050331 1.68% 
平均 相对 误差 14.12% 平均 相对 误差 0.96% 
序列 GM(1,1) 模 拟 ”相对 误差 高 阶 AGO 模型 模拟 值 相对 误差 
8 \ \ \ \ 
1 2.829706 182.97% 0.747971 25.20% 
6 3.247475 45.88% 6.107992 1.80% 
3 3.726922 24.23% 3.125645 4.19% 
4 4.277153 6.93% 4.009595 0.24% 
平均 相对 误差 65.00% 平均 相对 误差 7.86% 
这 里 同样 是 出 于 对 比 效果 的 目的 , 暂 不 考虑 原 序 列 的 模型 前 检验 而 直接 建 模 (实际 上 它们 的 
1-AGO 序列 均 能 通过 GM(1,1) 建 模 前 检验 ， 而 它们 本 身 不 能 通过 建 模 前 检验 ) . 


表 2 十 分 清晰 地 展现 了 高 阶 AGO 模型 的 强大 . 理论 上 来 说 , 正 是 因为 (90) 中 对 


2 一 Amx0(6) 的 引进 , 它 的 正 负 性 使 得 GM(1,1) 建 模 机 秆 
更 精准 地 还 原 了 序列 . 


4 序列 的 预测 


目前 为 止 , 我 们 运用 高 阶 AGO 生成 和 GM(1,1) 模 型 ， 


的 模拟 . 接 下 来 便 是 对 序列 发 展 的 预测 ， 正 如 传统 GM(1,1) 所 做 的 那样 . 


1 中 单调 性 的 缺陷 得 到 了 弥补 ， 从 而 


以 及 模拟 还 原 方法 实现 了 对 序列 


设 待 预测 项 为 rw (KE 十 D). 根据 (6)、(07) 式 , 待 预测 项 与 其 一 项 存在 关系 


在 (8) 中 x 了 (Kk 十 1])、 


Xm(k) 的 修正 项 ,也 是 已 知 量 


x(KE 十 D) 一 
二 (1 一 O44) x Kk 1 


x (FE) =xO (K+1) 
FD — (1—6) x (Kk) 


(8) 


x 中 (Kk) 是 已 知 量 ,是 通过 对 x 中 (1) 建 模 后 ,进行 IAGO 得 到 的 . 64 是 


.所 以 (8)F 


唯一 的 未 知 量 是 64 141. 


但 由 于 它 来 源 于 相对 误差 ， 


客观 上 并 不 存在 内 在 规律 , 所 以 不 建议 对 其 进行 建 模 预测 , 一 个 可 行 的 办 法 是 根据 6 的 历 


中 的 最 大 值 . 反之 , 妇 
影响 较 大 ,可 以 考虑 仅 对 预测 项 前 两 三 项 蔗 
进行 调整 ,按照 灰 数 理论 ! 
1,6,3,4} 


下 面 以 {8， 


并 行 高 


行 高 阶 AGO 模型 预测 示意 


史 信 息 对 其 进行 赋值 ， 这 也 是 灰 数 理论 中 常用 的 办 法 99 
如 果 希 望 保守 建 模 , 得 到 比较 小 的 预 
I 果 得 到 比较 大 的 预 


让 


测 值 ， 可 以 取 61+41 二 max《6:}, 也 就 是 654 历史 信息 
测 值 ， 则 可 以 取 6441= 
的 64 求 加 权 平 均 
赋值 的 技巧 结合 实际 情况 进行 赋值 . 
阐 列 为 例 进 


min{64}. 如 果 数 据 受 新 息 


下 赋值 . 总 之 ， 可 以 根据 预测 需求 


表 3 高 阶 AGO 模型 预测 示意 图 


序列 模拟 值 相对 误差 

8 \ \ 

1 0.747971 25.20% 

6 6.107992 1.80% 

3 3.125645 4.19% 

4 4.009595 0.24% 

平均 相对 误差 7.86% 
小 值 预测 2.071898 \ 
大 值 预测 9.739813 \ 
均值 预测 6.493338 \ 


表 中 均值 预测 是 指 取 5x 1 = avreage {64}. 
表 3 展示 了 高 阶 AGO 模型 的 三 种 预测 结果 . 在 对 原 序 列 的 模拟 上 , 平均 相对 误差 仅 有 
， 7.86%， 而 这 个 序列 本 身 无 法 通过 准 指数 律 检 验 从 而 建立 GM(1,1D) 模 型 . 从 表 2 可 以 看 到 , 直 
接 对 其 进行 建 模 的 话 , 平均 相对 误差 将 达到 65%. 


5 总 结 与 讨论 


a 


至 此 , 我 们 通过 对 振荡 序列 的 1-AGO 序列 建立 GM(1,1) 模 型 ， 再 运用 第 2 节 介 绍 的 模拟 
还 原 方法 还 原 到 对 振荡 序列 的 模拟 , 建立 了 基于 高 阶 AGO 生成 的 GM(1,1) 模 型 . 一 般 来 讲 ， 
由 于 多 次 IAGO 操作 会 造成 误差 的 累积 , 根据 文献 17 关 于 还 原 误差 的 讨论 , 设 传统 GM(1,1) 
模型 的 平均 误差 是 ce， 如 果 存 在 r 阶 AGO 模型 还 原 回 原 序列 的 方法 ， 那么 其 不 引入 修正 项 时 
的 平均 相对 误差 将 达到 2"e. 也 就 是 ， 随 着 建 模 阶 数 的 上 升 而 呈 2 的 指数 增长 . 但 是 在 引入 
了 修正 项 后 , 理论 上 会 将 这 个 误差 减 小 到 2"-1e. 虽然 在 > 很 大 的 时 候 这 个 误差 仍然 很 大 , 但 
振荡 序列 的 1-AGO 往往 已 经 能 够 通过 GM(1,1) 模 型 的 建 模 前 检验 ,实际 应 用 中 7 通常 取 2， 
特殊 情况 取 3, 理论 上 的 还 原 误 差 累 积 很 小 . 所 以 , 对 振荡 序列 建立 2 阶 AGO 生成 GM(1,1) 
模型 已 经 能 够 起 到 比较 好 的 效果 ,所 以 不 建议 对 原 序列 的 2 阶 以 上 的 AGO 序列 进行 建 模 ， 
除非 能 够 找到 这 种 情况 下 能 很 好 地 还 原 模拟 的 方法 ， 正 如 本 文 第 2 节 所 做 的 那样 . 

这 种 方法 相 较 于 利用 GM(1,1) 模 型 对 振荡 序列 建 模 目前 流行 的 三 种 方法 一 一 GM(1,1) 模 
型 群 、 残 差 GM(1,1) 模 型 和 区 间 预 测 而 言 ， 要 更 加 简便 和 深刻 . 所 谓 模型 群 , 就 是 将 振荡 序列 
分 成 许 许 多 多 的 小 区 间 , 然后 逐一 进行 建 模 08; 残 差 模型 本 质 上 是 一 种 二 次 建 模 , 即 进行 一 
次 建 模 后 ， 对 误差 序列 再 进行 一 次 建 模 的 方法 420; 区 间 预 测 是 对 区 间 的 极 大 值 和 极 小 值 分 
别 建 模 ， 从 而 能 够 找到 竺 预测 区 间 的 振荡 范围 ， 再 利用 灰 数 理论 中 的 白化 系数 ， 对 振荡 范围 
内 的 值 进行 赋 权 的 办 法 已 31. 而 本 文 介绍 的 高 阶 AGO 模型 ， 比 区 间 预 测 更 加 具有 准确 性 ， 
比 模型 群 的 办 法 更 加 方便 简洁 . 跟 残 差 模 型 原理 相似 , 但 在 序列 还 原 阶 段 就 对 误差 序列 进行 
了 处 理 , 无 需 再 次 建 模 这 样 繁琐 的 操作 . 
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The GM (1,1) Model Based on High Order AGO Generation 


Yuan Tu, Jiayu NI 
( College of Physics, Sichuan University, Chengdu 610065, China ) 


Abstract: In recent years, Grey System theory has been widely used in various fields. Among the Grey 
Prediction Models, the GM (1,1) model is the core and foundation. However, due to the exponential 
time response sequence, GM (1,1) model is difficult to simulate the oscillation sequence, and the 
oscillation sequence can not pass the GM (1,1) pre modeling test. These factors weaken the application 
of GM (1,1) model. The paper establishes GM (1,1) model for the 1-AGO sequence of oscillation 
sequence, by using the advantage of monotone sequence simulation of GM (1,1) model. Then, the 
IAGO operation with correction term is introduced to restore the simulation of oscillation sequence. 
Finally, an improved GM (1,1) model is established to make up for the defect of the traditional GM (1,1) 
model in the simulation of oscillation sequence. 
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